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We address a long standing issue and determine the 
decorrelation efficiency of the Hybrid Monte Carlo algorithm 
(HMC), for full QCD with Wilson fermions, with respect 
to vacuum topology. On the basis of five state-of-the art 
QCD vacuum field ensembles (with 3000 to 5000 trajectories 
each and ^--ratios in the regime > 0.56, for two sea quark 
flavours) we are able to establish, for the first time, that HMC 
provides sufficient tunneling between the different topological 
sectors of QCD. This will have an important bearing on the 
prospect to determine, by lattice techniques, the topological 
susceptibility of the vacuum, and topology sensitive quantities 
like the spin content of the proton, or the rj mass. 

PACS numbers: 11.15.Ha, 12.38. G, 83.20.JP 



I. INTRODUCTION 

In the past, a reliable access to topological quantities 
of the QCD vacuum by methods of Lattice Gauge Theory 
has turned out to be a real challenge to the entire lattice 
approach. 

Considerable progress has been achieved recently in the 
understanding of how to extract topological observables 
which are inherently continuum quantities from the dis- 
crete lattice. It was shown ]jj that the field theoretical 
0] and the geometrical [|| definitions of the topological 
charge yield — when suitably renormalized — equal values 
for the topological susceptibility, \i m the continuum. 
For practical purposes, however, we prefer the former 



definition, as it offers superior quality of the stochastic 
signal [0. 

This being settled, there remains the severe problem: 
any trustworthy stochastic sampling method must qual- 
ify to sufficiently decorrelate the members of the Monte 
Carlo time series with respect to the observables of in- 
terest; and the folklore is that the topological charge is 
the slowest of them all. 

This nuisance is particularly aggravated in presence of 
dynamical fermions in QCD as they induce long range 
interactions, in form of the determinant of the fermionic 
matrix M inside the path-integral. Standard Hybrid 
Monte Carlo algorithms (HMC) deal with this non-local 
problem by recourse to a stochastic Gaussian representa- 
tion of det(M) ||]. Such treatment implies the repeated 
costly solution of large systems of linear equations which 
puts narrow limits on the affordable sample sizes, even on 
Teracomputers. Thus autocorrelations on the HMC time 
series are an issue of prime importance in judging the 
statistical reliability of estimates from HMC in realistic 
settings. 

A failure of stochastic methods to scan the topologi- 
cal sectors of the QCD vacuum would be a real blow to 
lattice gauge theory, as topology is known to be of great 
importance to elementary particle physics: think of the 
Witten-Veneziano explanation of the rj' mass, the relation 
of the axial vector current divergence to the topological 
charge density through the Adler-Bell-Jackiw anomaly, 
or the role of instantons in the structure of the vacuum. 

For staggered fermions, indications have been given for 
quite insufficient tunneling rates of the topological charge 
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at very small mass parameters, both for two and four 
flavors |J|]. 

Some of us [§1 have performed simulations with four 
flavours of staggered fermions using the HMC. At (3 = 
5.35, on a 16 3 x 24 lattice, the time history exhibits little 
mobility, over 450 trajectories if the quark mass is cho- 
sen as ma — 0.01 corresponding to ^ = 0.57(2). Much 
shorter autocorrelation times, of the order of 10 trajecto- 
ries, are found at ma — 0.02 and ma = 0.05. These quark 
masses correspond to = 0.65(3) and & = 0.75(1). 

L rap v ' rap v ' 

For full QCD with dynamical Wilson fermions trust- 
worthy results for the topological decorrelation efficiency 
of fermion algorithms are missing in the literature. In 
this note we report on an investigation that will fill this 
gap. 

Our results are based on the high statistics samples 
as obtained by the SESAM |[o| and TxL [|l| projects 
with HMC on APE 100 hardware in Italy and Germany. 
The present sample consists of three time histories on 
lattices of size 16 3 x 32 @ at k = 0.156, 0.157, and 
0.1575, corresponding to three intermediate ratios 
of 0.839(4), 0.755(7), and 0.69(1), respectively. At the 
coupling [3 — 5.6, the scale computed from the p mass 
is a -1 = 2.33(6) GeV [^3[ , after chiral extrapolation. 
In the framework of the T%L project, SESAM's set has 
been complemented by two ensembles of configurations 
on 24 3 x 40 lattices, again at (3 = 5.6, with n = 0.1575 
and 0.158, the latter value corresponding to the lightest 
sea-quark mass attempted till now for Wilson fermions, 
with ^ = 0.56(2) II. 

So far, the SESAM sample has been partially analysed 
for integrated autocorrelation times of standard observ- 
ables only, like plaquette, extended Wilson loops and 
octet hadron masses jyj. The autocorrelation time, Tint 
of these observables appears to be bounded by the one 
of the lowest eigenvalue of the Wilson fermion matrix, 
the bound being of size 15 to 30 in units of molecular 
dynamics time. In this sense, we can start with sets of 
about 200 such 'decorrelated' configurations to search for 
tunneling through the topological sectors. 



II. TWO WAYS TO MONITOR FOR 
TOPOLOGICAL CHARGE 

The topological charge Q is the integral over space- 
time of the topological charge density, Q(x), which is 
related via the Ua(X) anomaly of the flavor singlet axial 
current, d^J^x) = 2N f Q(x). 

To monitor for the topological content along the time 
histories, we use two independent methods to estimate 
Q: 

1. The gluonic calculation, which starts from the 
charge density 



For the lattice definition of Q we follow Ref. 
Local fluctuations of the gauge fields are suppressed 
by the standard cooling procedure [Q] just as in 
Ref. §. 

2. The fermionic evaluation follows the Smit and Vink 
proposal J^] which is inspired by the Atiyah-Singer 
theorem: 



= mKp(Tr(7 5 G))j/. 



(2) 



Here G represents the quark propagator in the 
background gauge field U with quark mass m. up 
is a renormalization constant. No cooling is applied 
here. 

We comment that both methods are in principle prone 
to characteristic systematic errors: while the gluonic ap- 
proach might suffer from loss of instantons due to cooling, 
the fermionic monitor might be affected by 'dislocations'. 
It would therefore be gratifying to see them agree in their 
monitoring functionality. 

To our knowledge, the Smit- Vink proposal has never 
been really exploited. Computing the trace in Eq. (j^) re- 
quires the application of stochastic estimator techniques. 
Such estimates in general turn out to be extremely noisy. 
However, thanks to recent progress in noise reduction 
methods achieved by some of us [ |l5j , we are now in the 
position to extract rather accurate signals for Tr^G). 



III. RESULTS 

With the procedures decribed in the previous section 
we are sufficiently equipped to monitor for the decorre- 
lation of the topological charge in the HMC process. 

Fig. (l|) shows the time series for Q as measured on our 
ensembles. The first three graphs contain 200 configura- 
tions each. These ensembles were taken from contiguous 
series of trajectories of length 5000, produced by SESAM 
on 16 3 x 32 lattices. The configurations are separated by 
25 units in absolute molecular dynamics time. In the first 
three graphs, the results from the gluonic and fermionic 
measurements are super-imposed. We find them to agree 
nicely which demonstrates that the systematic errors of 
both methods are sufficiently under control. 

First inspection shows many tunneling events of the 
topological charge between the distinct sectors, for all 
three sea quark masses on the 16 3 x 32 lattices. There 
seems to be a tendency that the fluctuations decrease 
with increasing k. In an attempt to quantify these fea- 
tures we introduce the mobility of tunneling, Dj, as the 
average weighted number of tunneling events measured 
on configurations which are separated by a distance of d 
along the time history: 
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with Q being the integer value of the measured quantity 
Q and N the number of measurements. The resulting 
numbers are given in table ||. They reflect the tunneling 
rates and thus the decorrelation capability of HMC, when 
approaching the chiral limit. 

It should be noted, however, that Dd depends on the 
volume for two reasons: Dd has not been normalised to 
the amplitude of the fluctuations and our HMC parame- 
ters have been retuned when the system size was changed. 
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FIG. 1. Time series of topological charge. We plot series 
from 3 SESAM and 2 TxL ensembles. The last canvas is a 
quenched HMC simulation at (5 = 6.0. 



With the entire time histories of our HMC runs be- 
ing archived we are able to postprocess all trajectories 
from the SESAM and T%L samples. For k — 0.1575 — the 
smallest SESAM mass — we have computed the topolog- 
ical charge on a contiguous segment of > 5000 trajec- 
tories for every second trajectory, see Fig. (||). We find 
the topological charge to change frequently on that time 
scale. This is illustrated in the zoom of Fig. (|J), which 
shows considerable fluctuations on that time scale. 

The series for k = 0.1575 is long enough to allow for 
the computation of the autocorrelation function, as plot- 
ted in Fig. (||). We can apply an exponential fit to its 
'slowest' mode. The fit reveals the exponential autocor- 
relation time of r exp = 80(10). The integrated autocor- 
relation time from the lower diagram in Fig. (^|) turns 
out to be Ti n t = 54(4). These numbers should be com- 
pared to the upper bound of n n t as estimated from other 
observables, maxTi nt = 42(4) [Q. 

The tunneling behaviour is reflected in the histograms 
of Fig. (@). The errors quoted have been estimated by 
naive jackknifc. The distributions are well peaked at Q — 
0. With 200 entries, though, the symmetry is not yet so 
well established. 

Histogramming the fine scan of Fig. (j^), we arrive at 
the fourth histogram which comprises 2500 entries. We 
find a smooth Gaussian shape to emerge, see Fig. (0). 
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FIG. 2. Fine scan of of Q for re = 0.1575 on the 16 3 x 32 
lattice. The second canvas is a zoom. 

Finally, we turn to the volume dependency. On the 
24 3 x 40 system, at re = 0.1575, the time history looks 
rather much alike to the one on the 16 3 x 32 lattice, apart 
from an apparent increase in amplitude. The mobility at 
a separation of 25 units of molecular dynamics time, D25, 
turns out to be twice as large as for the small system, 
see table |. Naively, we would expect that the mobility 
scales approximately with the square root of the volume, 
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leading to a factor of two. However, we remark that on 
the large system, we had to re-adjust the length of the 
molecular dynamics trajectory to T — 0.5. Thus, for the 
chosen separation of T = 25, we refresh momenta and 
fcrmions twice as often as on the smaller system. 

Going more chiral, at k = 0.1580, our present sam- 
ple size is not sufficient to make definite statements, as 
the expected symmetry around (Q) = is not yet estab- 
lished. 

The last entry in Fig. (|l|) refers to a HMC series in 
the quenched case at j3 = 6.0, on a 16 4 lattice. The 
reduced volume amounts to a smaller amplitude. The 
series appears to be less ergodic on the time scale of 4000 
trajectories, compared to the three full-QCD runs on the 
16 3 x 32 lattice. The mobility in the quenched case is 
observed to take the value 1.2. 



IV. DISCUSSION AND OUTLOOK 

We find that HMC is capable to create sufficient tun- 
neling between the vacuum sectors when dealing with 
dynamical Wilson fermions for — > 0.56. 
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FIG. 3. Autocorrelation function and integrated autocor- 
relation time of Q for k = 0.1575 on the 16 3 x 32 system. 



This is seen with both, gluonic and fermionic determi- 
nations of Q which proves their reliability as monitoring 
devices. 
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FIG. 4. Histograms of the topological charge. 
TABLE I. Mobility D 25 , see eq. |. 
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The autocorrelation rates as determined from various 
non-topological quantities are in accord with the autocor- 
relation time of the topological charge. It is interesting 
to note that the high frequency mobility of HMC with re- 
spect to topology is important for achieving the expected 
symmetry of the topological charge distributions. 

When approaching the chiral regime, at k = 0.1580 or 
^ = 0.56, we would expect to need a sample size of 
order 10.000 to achieve sufficient ergodicity with respect 
to topology. This is definitely in the reach of teracom- 
puting. 

It is evident that the fermionic force plays an impor- 
tant role for the tunneling efficiency of HMC: in the 
quenched setting, HMC appears to be less efficient in 
driving the system through the topological sectors. In or- 
der to understand this effect in more detail further studies 
are needed. 

At this point, with the available sample, we are in the 
position to study hadronic properties related to topology, 
like the proton spin content and the topological suscep- 
tibility for —-values down to 0.69. Work along this line 
is in progress. 

Another line of research being followed by some of us 
fl6|| is the tunneling efficiency of alternative stochastic 
algorithms, such as the multibosonic method fl7j]. 
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